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FOREWARD 


The  dimensionless  curves  presented  in  this  report  are  intended  for 
use  in  routing  the  dambreak  flood  down  a  dry  prismatic  channel.  Three 
important  properties  of  the  flood  wave  may  be  determined:  Time  of 
arrival  of  the  wave  front*  maximum  flood  depth,  and  time  of  maximum  flood 
depth.  The  routing  curves  were  prepared  from  the  results  of  numerical 
simulation  experiments  which  solved  a  dimensionless  form  of  the  St.  Venant 
equations.  The  propagation  of  the  flood  wave  tip  along  the  dry  valley  floor 
was  calculated,  as  well  as  propogation  of  the  negative  wave  as  it  was 
reflected  from  the  upstream  boundary  of  the  reservoir.  It  was  possible 
to  express  the  results  in  the  form  of  dimensionless  graphs  because  valley 
geometry  was  expressed  as  a  simple,  prismatic  channel. 

In  utilizing  these  graphs,  it  is  necessary  to  transform  the  irregular 
natural  cross  sections  into  one  representative  prismatic  section.  The 
important  properties  which  must  be  preserved  are  storage  and  conveyance. 
Three  basic  forms  of  cross  section  shape  are  permitted:  rectangular, 
triangular  and  parabolic.  The  equation  for  expressing  the  properties 
of  each  of  these  forms,  as  well  as  other  equations  necessary  to  utilize 
the  dimensionless  graphs,  are  presented  in  Annex  I  at  the  end  of  this 
volume. 

This  report  culminates  a  line  of  research  that  started  with  the 
project  "Numerical  Techniques  for  Routing  the  Dam  Break  Flood  on  a  ury 
Bed"  In  which  Drs.  Theodor  Strelkoff  and  John  Sakkas  developed  a  computer 
program  for  calculating  the  advance  of  a  dam  break  flood  in  a  dry  prismatic 
channel.  In  a  subsequent  research  project,  "Flooding  from  Ruptured  Dams," 
the  same  investigators  extended  their  work  to  include  a  solution  of  the 
non-dimensional  form  of  the  St.  Venant  equation.  This  basic  work  was 
extended  to  produce  the  routing  curves  presented  in  this  report.  In  view 
of  the  many  uncertainties  involved  in  routing  the  dam  break  flood,  the 
fact  that  a  reach  of  a  river  has  to  be  transformed  into  a  prismatic 
section  is  not  an  overwhelming  simplification. 


ABSTRACT 


Dams  are  subject  to  failure  and  the  damages  produced  by  the  result¬ 
ing  flood  are  extensive.  Alleviation  or  prevention  of  the  damaging 
effects  requires  the  knowledge  of  the  flood  characteristics.  These 
include  primarily  the  time  of  arrival  of  the  flood-wave  front,  the 
maximum  flood  level  and  the  time  at  which  the  maximun  flood  level  occurs 
after  dam  failure. 

Using  a  rational  computation  technique,  based  on  the  method  of 
characteristics,  dimensionless  graphs  of  the  aforementioned  flood 
characteristics  were  prepared  for  a  prismatic  channel  of  general  para¬ 
bolic  cross  section  and  several  values  of  the  parameters  involved.  These 
values  were  selected  to  cover  the  practical  range  of  the  field  conditions 
affecting  the  magnitude  of  the  dam-break  flood.  Results  were  obtained 
for  such  periods  of  time  that  the  flood  peak  advanced  adeouately  down¬ 
stream  and  either  subsided  or  stabilized  considerably. 

KEY  WORDS:  Dams;  Floods;  Graphs;  Hydraulics;  Rivers;  Waves. 
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DIMENSIONLESS  GRAPHS  OF  FLOODS  FROM  RUPTURED  DAMS 


I.  INTRODUCTION 

Dams  constructed  to  form  large  water-storage  reservoirs  are  subject 
to  failure  for  several  reasons.  According  to  a  survey  (2)  foundation 
failure  and  spillway  inadequacy  count  for  about  two-thirds  of  all  dam 
failures.  Though  not  so  frequent,  failures  due  to  acts  of  war  and 
earthquakes  may  cause  more  serious  damage  since  they  usually  occur 
without  any  previous  warning. 

Experience  has  shown  that  floods  resulting  from  the  sudden  collapse 
of  a  dam  forming  a  large  water-storage  reservoir  are  catastrophic. 
Damages  are  anticipated  to  be  higher  in  the  future,  should  such  a  mis¬ 
hap  occurs.  This  will  be  due  to  both  larger  sizes  of  dams  being  built 
and  the  increase  in  industry  and  population  density.  Hunan  settlements 
and  the  industries  on  which  they  depend  usually  lay  in  the  flood  plain 
of  rivers  on  which  dams  have  been  constructed. 

Knowledge  of  potential  inundation  areas  can  lead  to  the  establish¬ 
ment  of  rational  real-estate  zoning  criteria  and  procedures  for  the 
emergency  evacuation  of  populated  areas  below  the  dam.  In  fact,  this 
is  the  intent  of  California  Senate  Bill  896,  passed  into  law  by  the 
California  Legislature  in  1972.  This  law  requires  the  preparation  of 
potential  flooding  maps  for  any  dam  in  California,  the  partial  or  total 
failure  of  which  would  result  in  death  or  personal  injury. 

Rational  and  accurate  methods  for  the  preparation  of  such  maps  are, 
as  yet,  unavailable  for  general  use.  Work  done  at  the  University  of 
California  at  Davis  (3,  4,  5)  led  to  the  development  of  a  rational  com¬ 
putation  procedure  for  predicting  the  flood  wave  resulting  from  the 
sudden  and  total  failure  of  a  dam  in  a  prismatic  river  valley. 

In  this  report  a  brief  outline  of  the  theoretical  aspects  of  the 
procedure  is  given.  For  details  the  interested  reader  should  consult 
the  references  cited  above.  The  parameters  involved  and  the  variables 
which  express  the  practical  significance  of  a  dam-break  flood  are  mainly 
stressed.  Finally,  dimensionless  graphs  of  these  variables  for  certain 
values,  over  the  practical  range,  of  the  parameters  are  given.  With 
these  graphs,  a  quick,  approximate  computation  of  the  most  important 
aspects  of  a  flood  wave  following  the  failure  of  a  dam  is  possible.  An 
illustrative  example  is  given. 


II.  THE  MATHEMATICAL  MODEL 
1.  The  Governing  Equations 

Unsteady  flow  encountered  in  flood-wave  movement  in  a  prismatic 
channel  is  described  by  the  Saint-Venant  equations  (2,  6) 

,  4 . 55  i£  .  5  4  .  0  (la) 

3x  ax  at 

+  V  4  +  i  4  =  g  (SQ-Sf )  Ob> 

at  ax  ax 


A  bar  over  a  variable  indicates  a  dimensional  quantity,  unbarred  variables 
fire  dimensionless.  In  Eqs.l:  A(y)  «  cros§-§ectional  area  of  flgwj 
V(x,t)  =  Q/A  =  average  velocity  of  flow;  Q(x,t)  =  discharge;  y(x,t)  = 
depth  of  flow;  B(y)  =_top  width_of  flow;  x  =  distance  from  the  dam, 
positive  downstream;  t  =  time;  g  =  ratio  of  weight  to  mass;  Sq  =  channel 
bottom  slope  and  Sf  =  friction  slope.  In  this  work  Sf  is  evaluated  using 
the  Manning  relation 


where  Cy  =  1.0  in  the  metric  system  and  Cy  =  1-486  in  the  British  system 
of  units;  n  =  Manning  roughness  coefficient;  R  =  A/P  =  hydraulic  radius 
and  P  =  wetted  perimeter. 

2.  Non-dimensionalization 


_  Taking  characteristic  quantities:  L  for  distance,  T-  for  time, 

Y0  for  depth  and  V0  for  velocity,  the  following  dimensionless  variables 
are  defined 


x  =  x/L0  ;  t  =  t/T0;  y  =  y/YQ  ;  V  =  V/V0 
A  =  A/Yq  ;  B  =  B/Y0  ;  R  =  R/Y0 


Introducing  these  variables  into  Eqs.  1  and  defining 
L0  =  •  Y0  =  L0  S0 


V 


o 


Cu  R  2/3  s  1/2 
_  K0  ^0 
n 


(3a,b) 


(3c, d) 


2 


the  Saint-Venant  equations  take  the  dimensionless  form 


A3V  +  V3£  +  i£=0 
B  ax  3x  at 


IV  +  y3V  + 
at  vax 


*0*  3X 


1  n  w2,ROx4/3, 
L1  -  V  (tt)  J 


(4a) 

(4b) 


where  R0  is  the  dimensionless  hydraulic  radius  corresponding  to  depth 
y  *  Y0  or  to  a  dimensionless  depth  of  unity. 

For  simplicity  Y0  is  taken  equal  to  the  maximum  value  of  water  depth 
behind  the  dam.  Then,  for  a  given  channel_geometry,  bottom  slope  and 
roughness,  Eqs.  3  yield  the  values  of  V0,  L0,  To  and  F0. 

Equations  4  are  applicable  to  channels  with  non-zero  bottom  slope. 

3.  Channel  Geometry 

In  this  work  a  prismatic  channel  is  assumed  with  breadth  B  related 
to  depth  y  by  a  formula  of  the  type 

B  -  C  f  (5) 

where  C  and  M  are  constants.  Non-dimensionalization  of  Eq.  5  yields 

B-Cy*1  (6) 

where 

c  =  c  ij'1  (7) 

If  B0  denotes  the^breadth  corresponding  to  depth  y  =  Y0,  it  is  easily 
seen  that  C  =  B0/Y0,  the  relative  top  width  of  the  channel  cross-section 
at  the  dam  site. 


4.  Solution  of  the  Equations 

For  details  in  solving  Eqs.  4  and  for  the  initial  conditions,  upstream 
and  downstream  boundary  conditions  used  in  the  solution,  reference  should 
be  made  to  (3,  4,  5). 

In  brief  Eqs.  4  are  transformed  into  characteristic  form  and  solved 
numerically  over  the  irregular  grid  formed  by  the  characteristic  lines 
in  the  x-t  plane  using  a  predictor-corrector  scheme. 

The  water  behind  the  dam  is  assumed  to  be  at  rest  prior  to  dam 
failure  and  the  downstream  channel  is  dry.  A  generalized  Ritter  solution 
(3,  4)  is  taken  to  represent  initial  conditions  along  a  forward  charac¬ 
teristic  shrunk  to  the  point  x  =  0,  t  =  0  in  the  x-t  plane. 


The  front  of  the  negative  wave  propagating  into  the  still  water  in 
the  reservoir  defines  the  upstream  boundary  of  the  flow  field  where  con¬ 
ditions  are  the  same  as  the  undisturbed  initial  conditions.  Upstream 
boundary  conditions  are  determined  in  this  fashion  until  about  the 
negative  wave  front  reaches  the  upstream  reservoir  end  at  time  tE  given 
by 

tE  =  2F0  /TST  (3) 

Shortly  before  that  time,  when  a  depth  smaller  than  a  prescribed  depth  ym 
is  obtained  at  the  moving  upstream  boundary,  a  fictitious  stream  of  the 
same  depth  is  introduced  at  the  ultimate  location  of  the  upstream  boundary 
and  assumed  existing  thereafter.  The  velocity  is  determined  using  the 
backward  characteristic  relation.  Typical  values  of  ym  are:  y  =  0.01 
for  f0  >  0.5  and  ym  =  0.03  for  F0  =  0.025.  m 

In  the  region  very  close  to  the  wave  tip  where  y-*D  and  dy/dx  -*■ 
the  formal  numerical  solution  is  ver7  costly  to  apply.  Instead,  the  water 
surface  profile  is  determined  analytically  using  a  simplified  form  of  Eq. 

4b  suggested  by  Whitham  (7)  and  based  on  physical  considerations  of  the 
tip  region.  Typical  values  of  y  up  to  which  the  formal  computation  proceeds 
are  in  the  order  of  0.04. 

Solution  of  Eqs.  4  yield  values  of  depth  and  velocity  at  the  nodes 
of  the  characteristics  grid  in  the  x-t  plane.  Data  pertaining  to  stage 
and  discharge  hydrographs  at  various  locations  along  the  channel  are 
obtained  through  linear  interpolation  between  node  values. 
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111.  DIMENSIONLESS  FLOOD  GRAPHS 
1 .  Dam-break  Flood  Variables 


For  practical  purposes  it  is  desirable  and  adequate  that  the  time  of 
arrival  of  the  wave  front  as  well  as  the  maximum  flood  level  and  the 
time  of  its  occurrence  after  dam  rupture,  be  known  for  any  given  location 
downstream  of  the  dam.  In  some  cases  and  for  preliminary  investigations 
a  rough  estimate  of  the  above  factors  is  sufficient.  Under  these  condi¬ 
tions  approximate  results  can  be  drawn  quickly  from  properly  prepared 
graphs  representing  solutions  to  a  series  of  simple  cases. 


A  simplification  of  the  river  valley  geometry  leads  to  a  prismatic 
channel  of  general  parabolic  cross-section  as  defined  by  Eq.  5.  In  this 
case  Eqs.  4  contain  three  parameters,  namely  F0,  M  and  C,  besides  the 
dependent  variables  y  and  V  and  the  independent  variables  x  and  t.  Thus 
the  variables  y  and  V  are  functions  of  five  variables,  that  is 


y  =  q  (x»  t,  Fq,  m,  c 

V  =  r  (x,  t,  Fo,  M,  C 


) 


For  given  values  of  x,  F0,  M  and  C,  Eqs.  9  and  10  represent  the  stage  and 
velocity  hydrographs,  respectively. 

From  the  stage  hydrograph  and  for  practical  purposes,  one  is  mainly 


interested  in  the  maximum  value  of  depth,  y^,  occurrinq  at  time  tM  such 
that  in  Eq.  9 

HHt=tM=°  OD 

Then 

yM  35  q  (x»  tM,  F0,  M,  C)  (12) 

Solving  Eq.  11  for  tM,  one  obtains 

t  =  s  (x,  F0,  M,  C)  (13) 

*4 

Introduction  of  Eq.  13  into  Eq.  12  yields 

yM  =  u  (x,  F0,  M,  C)  (14) 

The  speed  of  propagation  of  the  wave  front,  W,  is  equal  to  the  flow  velocity 
at  the  wave  front.  If  xw  and  tw  designate  the  wave-front  location  in  the 
x-t  plane,  then 

W  =  =  r  (xw,  tw,  F0,  M,  C)  (15) 
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Integration  of  Eq.  15  and  solution  for  t#  yields 


(16) 


tw  =  v  (x^,  F0,  M,  C) 

The  left-hand  side  of  each  of  Eqs.  13,  14  and  16,  constituting  the  practically 
important  dam-break  flood  variables,  can  be  plotted  versus  x  with  C  as  a 
parameter  in  one  sheet  of  paper  for  given  values  of  fg  and  M.  In  the 
triangular  cross-section,  defined  for  H  =  1,  C  is  no  longer  a  parameter 
because  the  hydraulic  radius  through  which  C  enters  Eqs.  4  is  a  linear 
function  of  y  and  thus  C  is  eliminated.  In  this  case  F0  may  be  used  as  a 
parameter  to  distinguish  curves  In  the  same  plot. 

2.  Range  of  Parameters 

For  the  preparation  of  dimensionless  graphs  of  the  dam-break  flood 
variables  actual  conditions  are  assumed  to  vary  between  the  following 
extremes : 


Condi ti on 


Minimiro  value  Maximun  value 


Water  depth  behind  dam,  YQ 
Channel  bottom  slope,  S0 
Manning  roughness  coefficient,  n 
Cross-section  exponent,  M 
Relative  top  width,_C  =  B0/Yn 

Froude  nunber  l0  =  V0//  -a 

y  *o 


10  m  (32.8  ft) 
0.0005 
0.015 
0 
1 

0.025 


200  m  (656  ft) 
0.010 
0.150 
1 

50 

5 


The  extreme  values  of  F0  are  determined  from  the  combinations  of  Y0,  S0, 
n,  M  and  C  leading  to  extreme  values  of  V0. 


3.  Data  Acquisition 

For  preparing  the  dimensionless  flood  graphs,  data  were  obtained  for 
the  following  values  of  the  parameters  fo»  M  and  C  in  the  above  listed 
range: 

f0  =  0.025  0.10  0.50  1  2  5 

M  =  0  0.50  1 

C  ■  1  2  5  10  50 


The  wave-front  location  at  any  time  is  given  directly  by  the  solution. 
The  value  of  y^  in  each  stage  hydrograph  is  determined  as  the  maximun  of 
discrete  values  of  y  obtained  in  the  course  of  the  solution.  The  latter 
are  close  enough  to  Insure  practically  insignificant  deviation  from  the 
true  maximun.  The  time  at  which  y^  occurs  is  the  value  of  tfl. 

The  solution  progressed  in  time  until  values  of  y^  and  t^  were 
obtained  for  a  cross-section  at  a  distance  at  least  m?g  below  the  dam. 
Minimum  value  of  m  is  1000.  If  X|w  designates  the  abscissa  of  a  flood 
peak,  then 


(17) 


max  t  xM  }  >_  m  YQ 


or 


max  {  xM  }  _>  m  S0  (18) 

For  given  values  of  F0,  M  and  C,  Eq.  18  must  be  satisfied  for  all  possible 
values  of  S0  in  the  respective  range.  Hence  Eq.  18  is  rewritten  as 

max  {  Xfo  }  >_  m  •  max  {  S0  }  (S0  _<  0.010)  (19) 

In  this  case  one  obtains  from  Eqs.  3c,  d 

s„  -  9  (-^)2  ^T-  (20) 

Cu  Ro 

For  any  fixed  value  of  Y0  maximum  value  of  Sp  results  when  n  attains  its 
maximun  permissible  value  of  0.150.  For  such  value  of  n,Y0  is  varied  over 
its  range  and  the  maximum  value  of  S0  is  found.  If  the  latter  is  less 
than  0.010,  it  is  inserted  into  Eq.  19.  Otherwise  max  {  S0  }*  0.010. 

4.  Arrangement  of  the  Graphs 

The  dimensionless  flood  graphs  are  arranged  in  three  sets,  one  for 
each  of  the  dam-break  flood  variables,  namely  t*,  yM  and  ty.  For  easier 
identification  and  use  each  set  is  included  in  a  separate  Appendix.  In 
Appendix  A  the  dimensionless  graphs  for  the  time  of  arrival  of  the  wave 
front,  t*,  are  given.  Figs.  A1  to  A15.  Appendix  B  contains  the  dimension¬ 
less  graphs  for  the  maximun  flood  level,  y^.  Figs.  B1  to  B15.  Finally 
in  Appendix  C  the  diemsnionless  graphs  for  the  time  of  occurrence  of  the 
maximum  flood  level,  tM,  are  found,  Figs.  Cl  to  Cl 5. 

In  each  Appendix  there  are  three  subsets,  one  for  each  value  of  the 
cross-section  exponent  M,  that  is  M  =  0,  0.50  and  1.  In  each  subset 
individual  sheets  are  arranged  in  ascending  values  of  F0. 

5.  Example  of  Using  the  Graphs 

An  approximately  prismatic  river  valley  is  400  ft  and  600  ft  wide 
at  distances,  respectively,  20  ft  and  86  ft  from  the  valley  floor. 

The  roughness  of  the  valley  is  estimated  at  n  =  0.070  and  the  slope  is 
approximately  S0  =  0.001.  Water  is  impounded  behind  a  dam  at  86  ft 
depth.  It  is  required  to  determine  the  flood  characteristics  in  case 
the  dam  instaneously  and  completely  collapses. 
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Solution.  Fitting  a  curve  of  the  form  of  Eq.  Si, through  the  boundary 
of  the  val  1  ey  one  finds  M  *  0.278  anjj  C  *  175  ft0*722  (dimensions  of  C 
are  ft  to  the  (1 =M)  power).  Taking  Y0  =  86  ft,  one  finds  by  Eq*  7  C  =  7 
and  by  Eqs.  3  :  L  =  86,000  ft  or  16.3  mi.,  V0  =  10.54  ft/sec,  T0  =  2.265 
hr  and  F0  =  0.20. 

Using  the  graphs  in  Figs.  A2,  A3  values  of  t*  at  the  desired  locations 
are  obtained  for  F0=  0.10  and  F0  s  0.50,  respectively,  and  K  *  0,  C  *  7. 

In  a  similar  way,  from  Figs.  A8,  A9  values  of  t*  are  found  for  F0  »  0.10 
and  F0  *  0.50,  respectively,  and  M  *  0.50,  C  =  7.  Linear  interpolation 
yields  the  values  of  t*  for  Fo  a  0.20  and  each  value  of  M.  Finally,  a 
second  interpolation  between  the  latter  yields  values  of  t*  for  F0  *  0.20 
and  M  =  0.278.  An  arrangement  of  data  in  table  form,  as  shown  in  Table  1, 
makes  the  task  easier. 

A  completely  analogous  procedure  using  the  appropriate  graphs  in 
Appendices  B  and  C  yields  the  values  of  y^  and  t^  as  shown  in  Tables  2 
and  3. 


TABLt  1 

Time  of  arrival  of  the  wave  front 


miles 

0.10 

0.50 

0.20 

0.10 

0.50 

0.20 

** 

ft mm 

6.00 

0.368 

0.25 

0.375 

0.281 

0.22 

0.28 

0.235 

0.255 

0.58 

8.15 

0.5 

0.40 

0.50 

0.425 

0.36 

0.40 

0.370 

0.395 

0.90 

16.30 

1.0 

1.02 

1.20 

1.065 

0.95 

1.00 

0.962 

1.008 

2.2P 

24.45 

1.5 

1.80 

1.90 

1.825 

1.65 

1.75 

1.675 

1.742 

3.94 

32.60 

2.0 

2.65 

2.85 

2.700 

2.43 

2.60 

2.472 

2.573 

5.83 

48.90 

3.0 

4.65 

4.875 

4.706 

4.12 

4.30 

4.165 

4.405 

9.98 

TABLE  2 

Maximum  flood  level 


■■hI{37TT7*9HH 

- PI - 

fr=TT5B 

h  = 

TT775-I 

- - 

X 

miles 

X 

*0 

fo 

nmmSi 

0.10 

0.50 

0.20 

0.10 

0.50 

0.20 

6.00 

0.368 

0.433 

0.395 

0.424 

0.524 

0.488 

0.515 

0.475 

40.8 

8.15 

0.5 

0.408 

0.378 

0.400 

0.500 

0.470 

0.492 

0.451 

38.8 

16.30 

1.0 

0.344 

0.325 

0.339 

0.440 

0.423 

0.436 

0.393 

33.8 

24.45 

1.5 

0.305 

0.295 

0.302 

0.405 

0.392 

0.402 

0.358 

30.8 

32.60 

2.0 

0.279 

0.269 

0.276 

0.378 

0.369 

0.376 

0.332 

28.6 

48.90 

3.0 

0.235 

0.230 

0.234 

0.340 

0.333 

0.338 

0.292 

25.1 

T 


TABLE  3 

Time  of  occurrence  of  maximum  flood  level 


)  Distance 

- jnnj - 

■■■■UaU&llflHI 

X 

miles 

n 

fo 

*0  = 

0.10 

0.20 

0.10 

0.50 

0.20 

tM 

tM,  hr 

6.00 

0.368 

0.58 

0.75 

0.622 

0.54 

0.70 

0.580 

0.599 

1.36 

8.15 

0.5 

0.78 

0.90 

0.810 

0.75 

0.90 

0.788 

0.798 

1.81 

16.30 

1.0 

1.60 

1.80 

1.650 

1.53 

1.70 

1.572 

1.607 

3.64 

24.45 

1.5 

2.45 

2.65 

2.500 

2.35 

2.50 

2.388 

2.438 

5.52 

32.60 

2.0 

3.40 

3.60 

3.450 

3.22 

3.40 

3.265 

3.347 

7.58 

48.90 

3.0 

5.45 

5.65 

5.500 

■■ 

5.05 

5.20 

5.088 

5.271 

11.93 
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APPENDIX  A 

GRAPHS  FOR  THE  TIME  OF  ARRIVAL  OF  THE  WAVE  FRONT 
(Figs.  A1  to  A15) 


TIME  OF 


DIMENSIONLESS  TIME 


4  6  8  10 


DIMENSIONLESS  DISTANCE  X/Lo 


17 


Fio.  A4 


DIMENSIONLESS  DISTANCE  X  /Lo 


TIME 


dimensionless  DISTANCE  x/L 


22 


Fig.  A9 


TIME  OF 

WAVE  FRONT  ARRIVAL 


*1 


28 


Fig.  AT 5 


APPENDIX  B 

GRAPHS  FOR  THE  MAXIMUM  FLOOD  LEVEL 
(Figs.  B1  to  B15) 


29 


MAXIMUM  FLOOD  DEPTH 


°A/"A  Hldaa  SS31N0ISN3WIQ 


DIMENSIONLESS  DISTANCE  x/L 


dimensionless  distance  x/L 


dimensionless  distance  s/L 


38 


dimensionless  distance  x/L 


dimensionless  distance  s/L 


dimensionless  distance  j/L 


DIMENSIONLESS  DISTANCE  x/L 


APPENDIX  C 
GRAPHS  FOR 

THE  TIME  OF  OCCURRENCE  OF  MAXIMUM  FLOOD  LEVEL 
(Figs.  Cl  to  CIS) 


45 


dimensionless  DISTANCE  x/L 


DIMENSIONLESS  TIME  TM/T0 


DIMENSIONLESS  DISTANCE  S/L0 


47 


Fin.  CP 


DIMENSIONLESS  DISTANCE  S/L 


DIMENSIONLESS  DISTANCE  3/L0 


DIMENSIONLESS  TIME  TM/T0 


2  4  6 

DIMENSIONLESS  DISTANCE 


l/L, 


50 


Fig.  C5 


2  4  6  8 


DIMENSIONLESS  DISTANCE  S/L0 


dimensionless  distance 


53 


Fig.  C8 


TIME  OF 

MAXIMUM  FLOOD  DEPTH 


57 


Fig.  Cl 2 
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Fig.  Cl 3 


ANNEX  I 


Suggested  Procedure  for  Developing  a  "Representative"  Prismatic  Channel 
1-1.  Introduction 

Routing  the  dam  break  flood  with  these  non-dimensional  oranhs  is  a 
three  step  process.  The  procedure  in  the  main  body  of  this  report  is  the 
final  step.  Tt  requires  two  previous  steps  as  follows: 

The  Irregular  natural  topogranhy  must  be  transformed  into 
a  "representative"  prismatic  channel 

and 

The  coefficients  that  are  required  to  express  the 
resulting  cross  section  shape,  flow  velocity,  and  flow 
depth  must  be  non-dimensional i zed. 

The  following  procedure  is  a  suggestion  for  accomplishing  these  two 
steps.  It  is  designed  to  preserve  the  potential  energy  slope,  the  head 
on  the  dam  at  failure,  the  reservoir  volume,  and  the  distance  from  the 
dam  to  the  point  of  interest.  Notation  is  the  same  as  that  used  in  the 
main  body  of  the  report. 

1-2.  Representative  Stream  Slope 

Plot  a  stream  bed  profile  from  the  unstream  end  of  th*  reservoir  to 
the  downstream  point  of  interest.  Fit  a  sinnle,  straight  line  through 
the  plotted  profile  In  such  a  manner  that  the  average  potential  enemy 
gradient  Is  maintained.  This  represents  SQ  as  illustrated  in  Pin.  T-l. 

1-3.  Representative  Reservoir  Length 

Mark  the  location  of  the  dan  boino  analyzed  on  the  stmam  bed  pro^'le 
and  plot  the  reservoir  elevation  for  conditions  at  the  tine  of  failure. 
Measure  the  head,  7  ,  from  the  pool  elevation  to  the  channel  bed  ei?Vatior. 
Plot  this  head  above  the  SQ  profile. 


above  mouth 


This 


Calculate  lQ  by  dividing  the  head,  YQ,  by  bottom  slope,  S^. 
distance  represents  the  length  of  the  reservoir. 


T-4,  Calculating  the  Representative  Prismatic  Cross  Sec  Hon 

By  Idealizing  the  shape  of  the  reservoir  as  approximating  that  of  a 
pyramid,  the  volume  may  be  calculated  v<1th  the  followinq  equation. 

Vol  =  (1/?)  •  *  •  L  (1-1) 

A  =  Cross  sectional  area  at  dam  for  depth  Y 
L  =  Reservoir  length  for  depth  Y  at  the  dam 
It  Is  necessary  to  express  cross  sectional  area  in  terms  of  the 
coefficients  utilized  In  the  main  body  of  this  report.  Starting  with 
equation  (5),  page  3,  and  integratino  leads  to  the  followinn 


(I-?) 


Since  the  L  in  equation  1-1  may  also  be  expressed  in  t^rms  of  Y,  the 
volume  equation  may  he  written  as 


_  C 

Vo1  "  5  .  (M+V)'  •'  S  ’  Y 
o 


(T-3) 


This  equation  has  one  independent  term,  Y,  and  one  dependent  term, 
Vol.  The  slope  parameter,  SQ,  is  known,  and  the  two  coefficients  C  and 
M  and  the  unknowns.  By  curve  fitting  through  ?.  points,  two  equations  may 
be  developed  as  follows. 


A 


at  point  1; 


Voli s  3  •  (W"r?ro  ‘  viH+Z  (I'4) 

at  point  2; 

Vo12 =  t“  "f4r)~  '«r  ‘  ^M+?  (T_,;) 

Solving  equations  (1-4)  and  (1-5)  yields 

M  =  (logfVoyVol^/logfY^))  -  2  (1-5) 

c  *  Vol,  •  3  •  (f»D  •  s/r9  0-7) 

Since  M  Is  dimensionless,  any  units  may  be  used  for  Vol  and  for  Y  In 
equation  (1-6).  Equation  (1-7)  requires  consistent  units  between  Vol  and 
Y,  however.  The  following  example  illustrates  the  use  of  equations  (1-5) 
and  (1-7). 

Actually,  any  two  arbitrary  points  may  be  selected  for  calculating  C 
and  M.  However,  In  keeplnq  with  the  objective  of  approximating  the  reservoir 
volume,  the  head  on  the  dam  at  the  time  of  failure  Is  selected  as  one  of  the 
points  and  g-  of  that  head  Is  used  for  the  other.  These  are  converted  to 
elevations  and  reservoir  storages  read  from  the  elevation  capacity  curve 
as  follows. 


Table 

1-1.  Data 

for  Calculating  Coefficients 

for  the 

Prismatic  Section 

Head 

Channel 

Pool 

Reservoir 

Volume 

Point 

(Ft) 

Elevation 

Elevation 

ac  ft 

10®  cu  ft 

1 

86 

1323 

1409 

26572 

1157.476 

2 

38 

1323 

1361 

4134 

180.077 

Substituting  these  data  into  equations (1-6)  and  (1-7) 
M  *  (Log(41 34/26572 )/Log( 35/86))  -  2.0 
M  *  0.278 

C  *  180.077(10®)  •  3  •  1.278  •  0.001/3O?,?7R 


C  -  173.92  ft 
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1-5.  Calculating  the  Non-Dimensional  Coefficients 


The  coefficient,  C,  is  non-dimensional  i7»vl  by  usino  the  head  on  the 


dam  at  the  time  of  failure,  Y  . 


c  =  c 


Y  M_1 
o 


f  I-R) 


=  173.9? 
C  =  6.99 


86 


9.278  -  1 


The  Froude  number,  FQ,  is  already  non-dimensional.  However,  it  must 
be  calculated  for  conditions  at  the  moment  of  failure.  This  requires  the 
following  equations. 

£  ■  v/TC 


,,  _  1.486  -  2/3  1/2 

n  Ro  5o 


(I-**) 

(1-10) 


Calculatina  the  hydraulic  radius  term  in  equation  ( I -1 0)  Is  difficult 
since  wetted  perimeter  Is  not  known.  The  following  interpolation  scheme 
is  proposed.  When  M  *  0,  the  section  is  rectanoular;  when  M  =  1 ,  it  is 
a  triangular  section.  Using  the  value  calculated  above  for  C,  calculate 
the  hydraulic  radii  for  M  =  9  and  for  M  *  1.  Then  Interpolate,  based  on 


the  true  value  of  M,  to  obtain  as  follows. 
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With  the  hydraulic  radius  now  determined,  the  dimensional  velocitv 
may  be  calculated  from  equation  (1-10). 

In  this  problem,  SQ  Is  equal  to  0.001 . 

VQ  «  •  602/3  •  0.001 1/2 


=  10.5  fps 

The  Froude  number  Is  calculated  with  equation  (1-9). 

f0  *  io.5/^r 

F l  =  0.20 

The  required  dimensionless  parameters  are  completely  determined  and 
summarized  below 


M  =  0.278 
Fo  =  0.20 
C  =  7 

The  length  term,  LQ,  used  to  nondlmenslonal Ize  distance  Is  calculated  as 

*  Vs„  <M1> 

L0  ■=  86/0.001  =  86,000  ft 

The  time  used  to  calculate  dimensional  time  from  the  nondimensional 
routing  curves  Is 

fo  =  KA 

=  86,000/10.5 
=  2.28 


Tables  1,  2,  and  3  In  the  main  body  of  the  report  show  the  use  of 
these  dimensional  ess  parameters  for  routing  a  dam  break  flood.  The 
calculated  depth  should  be  referred  to  the  actual  bed  profile  for  expressing 
results  In  terms  of  elevation. 
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